Opto-mechanical probes of resonances in amplifying microresonators 
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We investigate whether the force and torque exerted by hght pressure on an irregularly shaped dielectric 
resonator allow to detect resonant frequencies, delivering information complemental to the scattering cross 
section by mechanical means. The peak-to- valley ratio in the torque signal can be many times larger than in the 
scattering cross section, and, furthermore, depends on the structure of the resonance wave pattern. The far-field 
emission pattern of the associated quasi-bound states can be tested by the angular dependence of the mechanical 
probes at finite amplification rate. We relate the force and torque to the scattering matrix and present numerical 
results for an annularly shaped dielectric resonator. 
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I. INTRODUCTION 

Waves confined in iiTegularly shaped geometries of optical 
microresonators (such as micro-optical lasers made of semi- 
conductors micro-crystals 0]^ or laser dye droplets jsl) 
pose practical and theoretical challenges, as an intricate inter- 
ference pattern arises from the multiple coherent scattering off 
the confining boundaries. This is especially true at resonant 
conditions, when the multiple scattering results in system- 
atic constructive interference. The most direct probes of these 
microresonators are scattering experiments: The systems are 
illuminated with a coherent light source, and the scattering 
cross section is detected. The resonant peaks observed in the 
cross section are related to quasi-bound states (found at com- 
plex energies or frequencies), which can be observed as the 
working modes of micro-optical lasers. Irregular geometries 
are favored because they offer a rich mode structure and per- 
mit highly anisotropic modes with well-defined directed emis- 
sion 0]. The resonance pattern in the total scattering cross 
section does not reflect this richness in the mode structure — 
according to Breit-Wigner theory, the weight and width of a 
resonance is determined by the life-time of the mode; the wave 
pattern itself is of no concern fs']. In this paper, we show that 
complemental information about the quasibound states is con- 
tained in the directly accessible opto-mechanical response of 
the system, which actually does distinguish between modes of 
different degree of anisotropy. 

The opto-mechanical probes of the resonances and the as- 
sociated quasi-bound states that we discuss in this paper are 
the force and the torque exerted on the dielectric microres- 
onator by the light pressure of an illuminating beam. We 
focus on the practically most useful case of irregular but ef- 
fectively two-dimensional geometries, and also allow for a fi- 
nite amplification rate. Over the past decade, opto-mechanical 
tools based on light pressure have found various applications 
in the manipulation of microscopic objects, for which usually 
simple shapes have been assumed, and precision-detection 



of the acting forces have become commonplace For in- 
stance, dielectric objects of micrometer dimensions have been 
brought into rotation by light pressure, and the acting torque 
has been determined from the rotation of the object in a vis- 
cous medium 

(3 S S [13 m El El ■ Typical torques are of 
the order of 10~^^Nm for micrometer-sized objects at 500nm 
wave length and lOmW intensity. These dimensions and op- 
eration parameters are also typical for optical microresonators 
and microlasers. The rotation rate depends on the viscosity 
of the ambient medium and is of order of several Hertz (typ- 
ical forces are of the order of 10 — lOOpN). In some of these 
experiments torques as small as ~ lO^^^Nm proved to be suf- 
ficient to induce rotation. The rotation technique can be used 
not only to determine the torque, but also to determine the vis- 
cosity once the tor que has been obtained by independent (e.g., 
optical) means I14llla.ll6ll . Micro-electromechanical systems 
1 17] have not been applied in this context so far, but in prin- 
ciple are also sensitive enough to detect the opto-mechanical 
response. 

We suggest to use the opto-mechanical response for ana- 
lyzing the internal wave dynamics in complicated geometries. 
We demonstrate that the mechanical response contains infor- 
mation which is analogous to the scattering cross section and 
the delay time (the conventional probes for resonances), but 
is sensitive to other, complemental aspects of the quasi-bound 
states, such as their degree of anisotropy (which is reflected 
by the torque). Consequentially, the mechanical probes help 
to distinguish between different wave patterns. At finite am- 
plification rate within the medium, they contain information 
on the far-field emission pattern. Our general argumentation 
is supported in a practical setting by numerical computations 
for the annularly shaped dielectric disk shown in Fig.^ which 
displays a multifaceted set of wave patterns due to its non- 
integrable classical ray dynamics |18| We use two dif- 
ferent numerical procedures, the wave-matching method (see 
e.g. Ref. Il9||) and the boundary element method |20j]. 

The organization of this paper is as follows: In Section Hll 
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FIG. 1 : The annular resonator used in the numerical investigations of 
this paper is composed of two circles with radii R and R' = 0.6 7?, 
and eccentricity S = 0.22 R. The refractive index is n = 1 outside 
the resonator, n = 1.8 in the annular region between the circles, and 
n = 3.3 inside the interior circle. 



we briefly discuss resonances, quasi-bound states, and con- 
ventional probes of their detection (the Wigner delay time and 
the scattering cross section) in the framework of scattering 
theory. Sectionlffllcontains the principal results of this paper 
Section ITlIAl provides the kinematical relations for the force 
and torque in terms of the information provided by the scatter- 
ing matrix. In Section lUlBl we compare the conventional and 
mechanical probes for the annular resonator and show that the 
anisotropy of the wave pattern is tested by the torque. Round- 
ing up the considerations, in Section lTlI Cl we discuss how the 
far-field emission pattern can be inferred by means of ampli- 
fication. Our conclusions are collected in Section IIVI The 
Appendix contains details on the derivation of the kinematical 
relations. 



We choose a circular region A of radius TZ > R, containing 
the resonator, and decompose the wave function in the exterior 
of A in the usual Hankel function basis, 

^ = J2 {at^H^^\kr) + a'^^'^H^^Hkr)) e""^ (4) 



where r and are polar coordinates. (Here and in the follow- 
ing, all sums run from — oo to oo.) The angularly resolved 
radiation in the far field is given by 
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Below the laser threshold, the expansion coefficients a 
of the outgoing wave are related to their incoming counter- 
parts Cm"' by linear relations 

(out) _ V <? ,0^'") 



(6) 



where the coefficients Smm' form the scattering matrix. The 
scattering matrix fulfills the time-reversal symmetry 



S—n, — m{ 1) 



(7) 



and is unitary for real k and n. For complex values, unitarity 
is replaced by 



(8) 



II. RESONANCES AND QUASI-BOUND STATES IN 
SCATTERING THEORY 

A. Scattering approach 

In this paper we apply standard scattering theory to the ef- 
fectively two-dimensional systems in question. We separate 
the two polarizations of the electromagnetic field, with either 
the electric field 



Ez = Re[exp(— io^t)-)/)] 
or the magnetic field 

Bz = c^^ Re[exp(— jwt)?/;] 



(1) 



(2) 



polarized perpendicular to the plane. In both cases, the com- 
plex wave function ip fulfills the two-dimensional Helmholtz 
equation 



(A + n^fc^)V' = 0, 



(3) 



where n is the position-dependent refractive index, k — lo/c 
is the wave number, and c is the speed of light in vacuum. 
Amplification is modelled by a complex refractive index with 
Imn < 0. 



B. Poles and quasi-bound states 

Quasi-bound states are found at complex values kc of k that 
permit nontrivial solutions 

^(out) ^ jjj 

the case of no inci- 
dent radiation. 
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It follows that the values kc are the poles of the scattering 
matrix, which (for real n) all reside in the lower complex plane 
as a consequence of causality. 

For the annular resonator, complex values kc are shown in 
the bottom panel of Fig.|2j and some typical wave patterns are 
presented in the right panels of Fig.|3 For this system with a 
reflection symmetry about the x axis, the quasi-bound states 
occur with even and odd parity |21], and further can be di- 
vided into whispering-gallery modes localized at the interior 
interface (class Wint) or at the exterior interface (class Wcxt), 
and more extended modes that we group into two classes 
A (inside the interior circle, these modes are of whispering- 
gallery character) and C (the remaining modes). Close to 
resonance in the complex A; -plane (fc 
scattering matrix can be approximated by \2, 
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where dc.m = (— l)'"ac,m due to time-reversal symmetry and 
we accounted for the potentially quasi-degenerate partner of 
opposite parity (indicated by the superscripts). 

In general, the quasi-bound states found for complex kc and 
real n can be transported to real values of fc ~ Re kc by set- 
ting Imn ~ ImkcKen/ Kekc, corresponding to an active, 
amplifying medium close to threshold |23]. Above thresh- 
old, poles formally move into the upper complex plane, which 
physically indicates instability, and the linear relation between 
a^"^'> and 0^°"'^ breaks down. Yet, in homogeneously ampli- 
fying media the lasingmodes are well approximated by the 
cold-resonator modes I24l HEL E6ll . The far-field emission pat- 
tern of the laser. 
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(11) 



is then given by Eq. Q, evaluated with the quasi-bound state 
Qc that wins the mode competition (the first mode that be- 
comes unstable). 



C. Resonances and their conventional probes 

The values kc are the poles of 5* where this matrix is singu- 
lar, and are reflected by resonances of the system at real values 
fc = Re kc- Conventional probes for resonances are obtained 
from the total scattering cross section 



,(m)|2 



and the weighted delay time |2 
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FIG. 2: (Color) The top panels show four quantities (in arbitrary 
units) as a function of k that probe for resonances in the annular 
resonator: The Wigner delay time tw, the angle-of-incidence av- 
eraged scattering cross section ao, the angle-of-incidence averaged 
force in forward direction Fo, and the variance of the torque var A'^. 
The lower panel shows the complex resonance wave numbers kc of 
quasi-bound states, classified by their wave patterns as described in 
the text: Wint (o), Woxt (A), A (□), C (o). Open symbols indi- 
cate modes of even parity, full symbols are modes of odd parity. The 
complex resonances are also indicated by the spikes in the top pan- 
els, located at fc = Re kc with height proportional to the life time 
1 / (—2c Im kc). Additionally, the resonances of Fig.|3ja-e) are indi- 
cated by the arrows at the very top. 



Clearly, the quantities a and t depend on the incoming wave, 
which we now specify as a plane wave coming from direction 
00 I2JI, corresponding to 



,(in) 



1 



-im{Tr/2+(po) 



(14) 



Both (7 and r then provide angularly resolved information as 
a function of ^q. A global characterization of the system is 
obtained by an average over the incident radiation direction 
00, giving 



CTO 
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(1 + El^" 



dS 



2 Re 5™™), (15) 
(16) 



where tw is known as the Wigner delay time. 

Both the delay time (of predominantly theoretical virtue) 
and the (more practical) scattering cross section display peaks 
at resonance, as is illustrated for the annular resonator in the 
two topmost graphs of Fig. |2] Note that the signal cto is, in 



general, of low contrast and, apart from a relatively strong 
background modulation, rather featureless. Still, the enhanced 
scattering of the light field at resonance promises a marked 
mechanical response, which we investigate in the remainder 
of this paper Indeed, the signal of the mechanical response 
will display a much better contrast for a set of systematically 
selected resonances. 



III. MECHANICAL DETECTION OF RESONANCES 
A. Kinematics 

In this subsection we provide general kinematic relations 
for a two-dimensional resonator in a light field, which in the 
following Subsections IIII BIIII CI will be used to characterize 
resonances and wave patterns. 

The mechanical forces exerted by the light field on a di- 
electric medium originate from the refraction and diffraction 
at the dielectric interfaces — ultimately, from the deflection 
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(and creation, at finite amplification) of the photons. The kine- 
matics in the combined system of Hght field and medium can 
be obtained from the conservation laws of total angular and 
linear momentum, which equate the torque and force acting 
on the medium to the deficit of the angular and linear mo- 
menta carried by the electromagnetic field into and out of the 
circular region A. (The center of this region is identified with 
the point of reference for the torque, and in our example is 
taken as the center of the exterior circle.) Note that the con- 
servation laws also hold in an amplifying medium, due to the 
recoil of each created photon. The kinematic relations hence 
follow from integrals of Maxwell's stress tensor |29| over the 
boundary of A. After some algebra (for details see the Ap- 
pendix), we find the time-averaged force and torque (per unit 
of thickness of the resonator) in the compact form ISflll 

J^x^l-^y — ^ /_^\"'m "m+1 "m "m+l J-iy^'' 
m 

^ = ^-^T.^{\4T^\'-\-irr)- (18) 

■m 

In absence of amplification and for plane-wave illumination, 
it is easily seen that the direction-averaged force and torque 
vanish from the unitarity constraints of the scattering matrix: 

{Fx+iFy)^^ — 2^ ^ •S'mm'S'm+l.m' = 0: (19) 

m,m' 




= 0. (20) 



Two simple quantities that do not vanish are the mean of the 
force component F|| — F^ cos 0o + Py sin 4>o in forward di- 
rection, 

(21) 

and the variance var = {N^)<j>o °f torque, 

^2 ' 

varA''=^^ mm' S^n_^Smn2Sm'n3^m'ni- (22) 

Here, the prime at the sum enforces the restriction ni + 713 = 

B. Resonances 

Because scattering is enhanced at resonance, we expect 
that var A^ and Fq are global characteristics of the resonances 
comparable to ctq and tw, while A^ and F|| provide angularly 
resolved information as a function of the incident radiation di- 
rection (j)o, analogously to a and r. The results in Figs.|2]and 
|3]demonstrate this promise to hold true. 

Figure 12] shows the global characteristics (Tq, t^y, va,rN, 
and Ff) for the annular resonator as a function of k. The plot 
demonstrates a clear correspondence of the resonant peaks in 



all four quantities. Evidently, each quantity probes another 
aspect of the resonances, such that their relative weights are 
different. As usual, the delay time tw displays the largest 
peaks for the very narrow resonances associated to long-living 
quasi-bound states. Presently, the longest-living states are 
those of class Wint, followed by those of class Woxt and A, 
while the states of class C are hardly visible here. The scatter- 
ing cross section ctq displays smaller peaks (sometimes, dips) 
on a modulated background, and does not provide a distinctive 
discrimination between the different states. 

The force and torque are sensitive to the wave pattern itself, 
and for the annular resonator display a marked response es- 
pecially for the quasi-bound states of class Woxt and A (for 
small k, also for class C). The peak-to-valley ratio in the 
variance of the torque (e.g., « 5 for the peak at /c « 6.2) 
can take much larger values than in the scattering cross sec- 
tion (where the ratio rarely exceeds 1 .2) for resonances with a 
rather anisotropic internal wave pattern. 

This principal conclusion of the present paper is supported 
in more detail by the five examples in Fig. |3] The wave pat- 
terns shown in the right panels are inhomogeneous to various 
extent, even though the inhomogeneity does not automatically 
translate into very anisotropic far-field emission patterns for 
these rather low values of Re kcR. Let us inspect the associ- 
ated resonance peaks in Fig. |2] For the quasibound state of 
Fig-EJa), force, torque, and scattering cross section detect the 
resonance with comparable contrast; the best signal is given 
by the delay time, which, however, must be remembered to 
be an inconvenient tool for practical considerations. Figure 
|3lb) pertains to a comparatively short-living mode (with large 
jlm/ccl). This mode is hard to detect by the conventional 
probes, but gives a very clear signal in the mechanical probes, 
especially, in the torque. For the case of Fig.|3jc), the best sig- 
nals are provided by the force and by the delay time, followed 
by the scattering cross section, while the torque is insensitive 
to this resonance. The case of Fig.|3ld) illustrates that reso- 
nances with both a short life time and a rather homogeneous 
wave pattern are hardly detected by any of the four probes — 
such modes, however, are also of subordinate interest in the 
practical applications. The mode of Fig. (Sje) gives a torque 
signal with a much higher contrast than the force or the scat- 
tering cross section; only the delay time signals the resonance 
to a similar extent. 

The most remarkable difference in Fig.|3]is between panels 
(b) and (d) — both modes of class C displayed there have 
a short life time, which inhibits their detection by conven- 
tional means, but the mode of panel (b) is detectable by opto- 
mechanical means, in keeping with the larger anisotropy dis- 
played by its wave pattern (note the regions of concentrated 
intensity above and below the interior circle). 



C. Far field 

In the practical application of micro-optical lasers, the 
anisotropy of the internal wave pattern is a desired design 
goal since it is the pre-requisite for directed emission in the 
far field. To which extent can the far field be inferred from 
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FIG. 3: (Color) The panels (a-e) show results for various com- 
binations of k and Imn tuned very close to resonance with the 
quasi-bound states of odd parity shown on the right, with (a) kc — 
6.009 - i 0.008 (class Wint), (b) fee = 6.228 - i 0.116 (class C), (c) 
kc = 7.554 - i 0.063 (class Woxt), (d) kc = 7.847 - i 0.145 (class 
C), (e) kc = 9.026-i 0.022 (class A). In the top graph of each panel, 
the angular dependence of the far field Ic (red) on the radiation di- 
rection (f} is compared to the angular dependence on the illumination 
direction </>o (note Ref. 1 28]) of the scattering cross section a (green), 
the weighted delay time r (blue), and the x-component of the force. 
Fx (purple) — the four lines are almost indistinguishable. In the bot- 
tom graph of each panel, the (squared) torque A*'^ (red) is compared 
to the interference term I'J^ (green), defined in Eq. <25> . The blue 
curve is N'^ evaluated at real n. All quantities are in arbitrary units. 



the angularly resolved information contained in N, F, a, and 
r? This rather innocent question turns out to be surprisingly 
subtle in view of the following observation: In absence of am- 
plification (real k and n), the angle-of-incidence averaged far 
field 

Jo(<^) = {im^o (23) 

is independent of cj). Here, we average over the angle of in- 
cidence (po since the quasi-bound states are defined without 
reference to any excitation mechanism. The angle of inci- 
dence (po enters the far field (|5} via Eq. (I14> . and after tak- 
ing the average the statement above follows from the unitarity 
of the scattering matrix. In other words, in the absence of 
amplification, the angle-of-incidence averaged far-field radia- 
tion pattern does not carry any intrinsic information about the 
quasi-bound states, even at resonant conditions. 

Can the far field of the quasi-bound states be inferred from a 
more sophisticated analysis of I{(p; (po), i.e., by taking the de- 
pendence of the angle of incidence into account? This would 
require the delicate task to discard the component of the light 
that is directly reflected at the first encounter of the bound- 
ary from the outside. The direct contribution is essentially 
independent of the quasi-bound states: the latter are deter- 
mined by a constructive-interference condition for reflection 
from the inside of the system while the directly reflected radi- 
ation never ever enters the medium. For instance, in our model 
system, the directly reflected wave component contains no in- 
formation on the interior circle, which, however, is crucial in 
the formation of all quasi-bound states. 

The situation changes significantly at finite amplification: 
Then, already the angle-of-incidence averaged intensity Io{(f>) 
is modulated, and is influenced by the quasi-bound states clos- 
est in k. At exact resonance in the complex fc-plane, approxi- 
mation (I10> of the scattering matrix entails with 
a large proportionality constant due to the resonant denomi- 
nator. Hence the contribution a''") of the incident radiation 
can be neglected, and I{4>) oc Idfp) defined in Eq. 
independent of the mode of excitation. Furthermore 

a(0o) cx r(0o) « \F,{^o)\ « 4(0 = 0o), (24) 

since again a*^""*' cx ac dominates over a^'"' in Eqs. il2\ . 
il3\ . (I17> . Equation (I24> entails a duality between the illumi- 
nation direction (j>o and the radiation direction cp. This duality 
relies on the time-reversal symmetry |28], which is incorpo- 
rated in Eq. dlOl l by the relation between Uc and ac. For a 
representative set of resonances at real k and complex n, the 
far field Ic{<P) is shown in the left panels of Fig.|3j along with 
the angular dependence of a, r, \Fx\, and N on the direction 
of incident radiation, (po. The wave pattern of the correspond- 
ing quasi-bound states is shown in the right panels. 

Due to the reflection symmetry of the annular resonator 
about the x axis, \Fy\ ^ \Fx\ while F\\ « cos{(po)Fx. 
The reflection symmetry also suppresses N compared to 
non-symmetric systems: The otherwise dominant contribution 
from vanishes for each given quasi-bound state, 

Sm'^l°'c,mP — 0, because of pairwise cancellation of the 
terms with opposite to. However, N is still enhanced com- 
pared to the non-resonant situation: A non-vanishing result 
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(of order \k — kc\~^) is obtained from the interference be 
tween the resonant state (with large coefficients) and the non 
resonant states (with moderate coefficients). In the typical 

{even 



case of quasi-degeneracy, Oc 
from Eq. ilO\ we obtain 



interferes with a^c"^'^\ and 



(25) 



This relation is indeed obeyed to good extent in the numerical 
computations (see Fig.|3}- Even at real n, N"^ (the blue curve 
in the bottom graph of each panel) roughly corresponds to 
I'^. In non-symmetric geometries the sum m|ac,mp ^ 
for a given quasi-bound state, and hence the proportional- 
ity N ^ Ic (and Fy ^ Ic) is restored (moreover, quasi- 
degeneracies are then lifted, and the system more easily is 
tuned to resonance with individual quasi-bound states). 



IV. CONCLUSIONS 

In summary, the force and torque exerted by light pres- 
sure on a dielectric resonator allow to detect resonances and 
help to characterize the wave patterns of the associated quasi- 
bound states. For anisotropic wave patterns that support a 
high angular-momentum transfer, the peak-to-valley ratio in 
the mechanical probes (notably the torque) exceeds by far 
the moderate values observed in the scattering cross section, 
which is notoriously insensitive to the wave pattern. 

We put our work in the context of directed transmission 
from micro-optical lasers, and took amplification into ac- 
count. Enhanced sensitivity to internal structure is also de- 
sired in several applications involving simpler passive or ab- 
sorbing media, such as cells with organelles or liquid drops 
polluted with inclusions. 

Since resonance provides a very effective scattering mech- 
anism by constructive interference, the typical forces and 
torques estimated for common micrometer sized dielectric 
resonators under typical radiation conditions are in the range 
of the opto-mechanical experiments mentioned in the intro- 
duction. In the viscous rotation experiments, the resonances 
will depend on the refractive index of the surrounding liq- 
uid, which may not be desirable for a precise characterization 
of the resonator, but also introduces an additional potentially 
useful control parameter. Detection by microelectromechani- 
cal systems offers the advantage to control the resonator ori- 
entation with respect to the incoming radiation. 

The numerical part of this work concentrated on the wave- 
optical regime, in which the wave length is not much smaller 
than the geometric features of the system and interference pat- 
terns are most complex. Some micro-optical lasers operate at 
smaller wave lengths (larger values of the wave number k) 
than accessed in our numerics. In this regime, semiclassi- 
cal relations between the internal wave pattern and the far- 
field emission pattern can be formulated, which also relate 
the directed emission desired for micro-optical laser to the 
underlying anisotropy of the internal wave pattern. The re- 
sults of this work hold the promise that in this semiclassical 
regime, modes with directed emission can be identified by the 



torque. The selectivity of the torque for anisotropic modes 
should be even enhanced for these larger values of k by the 
following mechanism: Isotropic modes frequently arise from 
a large collection of unstable ray trajectories and can be 
described as a superposition of random waves [32,1. As k is 
increased, more and more random-wave components become 
available, and hence the mechanical response is suppressed by 
self-averaging. Anisotropic modes are guided by just a few 
trajectories such that the self-averaging mechanism does not 
apply to them, and consequently they remain well detectable 
by their opto-mechanical response. 



APPENDIX: DERIVATION OF THE KINEMATICAL 
RELATIONS 

The calculation of the force 



F== / dSnT 
Is 



(A.1) 



and the torque 



N = [ dS (nT) X r (A.2) 
Js 



(where n is the unit vector in normal direction to the surface 
S) starts with the stress tensor 



1 



1 



1 



T ^ £o(E ® E - -EH) + — (B ® B - -BH). (A.3) 
2 Mo 2 

In the two-dimensional case, we integrate over the circle dA 
of radius TZ (the physical force and torque are obtained by 
a multiplication with the thickness of the sample), and it is 
natural to work in polar coordinates. 

Depending on the polarization, we insert the electromag- 
netic field 

B = Imw-le-''^*(7^-la0Ver - S^V^e^), J ^ 

B = c-iRee— *Ve., 1 
E = -cImw-le-*"*(7^-l50Ver - dri^e^), J ^ 

(A.4) 

where the wave function fulfills the Helmholtz equation (|3jl 
and is decomposed in the basis of Hankel function, Eq. @. 

After a time-average, the radial component of the stress ten- 
sor then is given in terms of 4' by 



(A.5) 



with Z = kU, dz dkr\r=n, k = uj/c, c = I/^moEo- 
(Note that k is real and rt = 1 on dA, since refraction and 
amplification is restricted to the resonator.) 

The force components can now be expressed as 

+zFy ^ '-^ j^J d<Pe^'^ [^-^m? - \dzi>? 

-IV'I' - ^<9^V'5zV'* - ^d^rdzi^ .(A.6) 
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Here we insert Eq. Q [indices m,j for V', indices m',j' for 
where the second index = (1) = out or (2) = in] 
and integrate over cj). Next, we use the identities 



dzHL 



^ rii -'-'m+li 



(A.7) 
(A.8) 



(here and in the following, we suppress the argument Z of the 
Hankel functions). This gives 



TTEoTZ 



ni{'m + 1) 



Z2 

m + 1 /TO 



Z \Z 



^^m^^m+l -'-'■(Ti-'-'m+l 
11 -'-'m+1 i -'-'m+1 



Between the square brackets, most terms cancel, giving 



(A.9) 



m,j,j' 

I -'-'m+l-'-'m -'-'m-'-'m+l 



(A. 10) 



With = and the identities 



4i 



An = -A22 - 2zImiJ^^|ii/,(2) = (A.ll) 



yli2 = A21 = 0, we arrive at the final expression (I17> . 
The calculation for the torque is less involved. We find 



N = 



; Bz im > TO 

k " ^-^ 



m— — 00 



a(i)*iy(2) , a(2)*^(i)A 

(2)^(2)/ , (1)^(1). 



(A. 12) 



The final result ([T8|l follows from = i/^^) and 



Imi7(?)(Z)i/«'(Z) = -|. 



(A. 13) 
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